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I. INTRODUCTION 



Currently the only missing piece of the Standard Model (SM) is the Higgs boson, which 
plays a vital role in this theory as it is assumed to be responsible for the mechanism of 
electroweak symmetry breaking (EWSB). Although only a single neutral CP-even Higgs 
boson is present in the SM, more than one Higgs boson may exist in several new physics 
models. For example, the two-Higgs doublet model (THDM), one of the simplest extensions 
of the SM, predicts the existence of three neutral and two charged Higgs bosons. In the 
CP-conserving THDM, two of the neutral Higgs bosons are CP-even and the remaining one 
is CP-odd. If CP violation appears in the scalar potential, the three neutral Higgs bosons 
may be a mixture of the CP eigenstates. 

The hopes to observe a Higgs boson at the CERN LEP collider faded away. The lower 
bound on the SM Higgs boson (H) mass m H > 115 GeV was obtained from studying the 
process e~e + — > ZH after combining the data taken by the four LEP experiments. The 
main Higgs boson production mechanism for a Higgs boson with mass around a few hundred 
GeV at a hadron collider, such as the Fermilab Tevatron and the CERN Large Hadron 
Collider (LHC), is through gluon fusion process. For ran < 2m#, the Higgs boson will 
decay predominantly into a bb pair, but detecting Higgs boson via this channel is challenging 
because of the large QCD backgrounds. In this mass range the rare decay modes H — > 77 
and H — > WW*, ZZ* are the most promising channels for detecting a SM Higgs boson. 
With enough integrated luminosity, the Tevatron Run-II is expected to offer the possibility 
of finding a Higgs boson with a mass up to 130 GeV through the associated WH or ZH 
production processes. At the LHC, the gluon fusion production mechanism will allow us to 
search for a Higgs boson via the H — > 77 decay mode for m H lower than about 130 GeV, 
whereas for m# ranging from 130 GeV to 180 GeV the H — > WW*, ZZ* decay modes could 
yield an observable signal jjj. 

Once a Higgs boson is detected, the next task will be to find out whether it is the SM 
Higgs boson or some other scalar bosons predicted by other theories. To this aim it is crucial 
to determine the properties of this particle, which may shed light on the EWSB mechanism. 
Therefore, apart from a precise determination of the Higgs boson mass (or masses if there is 
more than one Higgs boson), one of the main tasks of future colliders will be to determine 
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other properties of Higgs boson, such as its decay rates, its couplings to other particles, 
its properties of transformation under the discrete symmetries C (charge conjugation) and 
P (parity). Several methods have beenproposed in the literature to determine the CP 



properties of a Higgs boson at hadron 
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and photon |22j, |23j, 12J, l2Ia| colliders. Some of these methods 
rely on the study of certain energy and angular correlations of the decay products of the 
Higgs boson, whereas others concentrate on the change in its production rate due to the 
structure of its couplings to gauge bosons or heavy fermions. For instance, it was shown 
that polarized photon beams could be used to determine the CP property of a Higgs boson 
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of the Higgs boson were analyzed for the e e + — > ZH process [8]. In this work we show that 
the e~e + —> e~e + 4> process through VV (with V — Z, 7) fusion at the future International 
Linear Collider (ILC) could help us to determine the CP nature of the Higgs boson. It has 
long been known that for large center-of-mass energies (greater than about 500 GeV) the 

n n 

ZZ fusion channel becomes dominant in the e e + — > e e + (p process [26]. In Ref. (2JJ, it was 
shown how to measure the degree of CP violation in the Higgs sector via ZZ fusion process 
at a \/~S = 500 GeV ILC. As a complementary study, we will not consider CP violating cases 
in this work. We shall assume CP invariant interactions and study how to determine the 
CP property of the Higgs boson produced at the ILC with a higher center-of-mass energy 
(V^ = ITeV). Since at this energy, ZZ fusion contributes with 95% of the total cross 
section, as predicted by the SM, we choose to focus on the VV fusion channel alone. 

Rather than focusing on some specific model and calculating the respective loop contri- 
butions to the (j)VV vertex, our study is based on the general Lorentz invariant scalar-vector 
boson vertex. In general, three possible couplings (pZZ, (f>Zj and ^77, may give rise to ZZ, 
Z^i and 77 fusion processes, which at high energies will become the main source of Higgs 
boson production. As we shall see later, we found that the effects of the ^77 couplings are 
more effective for our task, as they can give large contributions to e~e + cj) production due to 
collinear enhancement. 

The rest of this work is organized as follows. In Section II we present a brief description of 
the general couplings. Section III is devoted to a Monte Carlo analysis for the e~e + —> e~e + (p 
process. Both polarized and unpolarized beams are considered and special attention is paid 
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to those kinematical distributions which are sensitive to both the strength of the anomalous 
(j)VV coupling and the CP property of the Higgs boson. Our conclusions will be presented 
in Section IV. Finally, the helicity amplitudes for the e~e + — > e~e + <fi process are presented 
in Appendix A for completeness. 



II. ANOMALOUS <pVV COUPLING 

The general Higgs-gauge boson vertex describing 077, fiZj and (f)ZZ couplings can be 
written as: 

V?w = i*+zzgr + z^f 1 (^-^Prft) + l ^e^p lp p 2a , (1) 

where p± and p v 2 are the momenta of the gauge bosons. The a^zz fi^term represents the 
dimension-4 coupling that appears in models such as the SM, the THDM and the MSSM. 
For simplicity, we will use the SM value a^zz = 2 M z /v in our study. 

The anomalous (pVV (with V = 7, V) couplings with dimension more than four could 
arise at the one-loop level in a renormalizable theory. Also, from the standpoint of the 



effective Lagrangian approach (ELA) [27, 28], these couplings can be induced by dimension 
6 (and higher) operators [3, For example, an effective Lagrangian using a SU{2) x U(l) 
scalar doublet that couples with the gauge bosons through dimension-6 operators would be: 



£37 = ^£(/«°«+/^) » ( 2 ) 

i 

where the CP-conserving operators Oj are given by (iul 
Obw = ^B^W^, 

o ww = ^w^w^s, 

O bb = $t B ^B^$, 

in which $ is the Higgs doublet, = ig' and = ig/2a a W^, with B^ u and W° the 
strength tensors of the gauge fields; g' and g are U(l) and SU(2) gauge couplings, respectively. 
As for the CP-violating operators O i} they are obtained from the above ones after replacing 
either of the two strength tensors by their respective dual |2i 



The simultaneous presence 



of both g$yy and g^vv for one boson would signal CP violation. As mentioned, we will not 
consider that case in this work, for this we refer the reader to Ref. j^jj. In our study we 
assume a SM-like CP-even Higgs boson h with anomalous couplings guvv and/or a CP-odd 
Higgs boson A of same mass but with couplings cjavv 1 We will also refer to the SM Higgs 
boson H for the purpose of comparison. Throughout this paper, we refer the production of 
the SM Higgs boson H as the one generated from the Born level dimension-4 HZZ coupling. 
Even though higher dimension anomalous operators may arise from loop corrections, we will 
assume their effects to be negligible. In contrast, for the production of the SM-like Higgs 
boson h, we include both the Born level dimension-4 HZZ coupling and the higher dimension 
anomalous ghvv couplings in our calculations. For the production of the CP-odd Higgs boson 
A, only the higher dimension anomalous cjavv couplings are included. 

A common approach for determining the CP property of Higgs boson is to examine the 
Z(p associated production (e~e + — > Zcj) — > f~ f + (fi) and analyze the angular distribution 
of the Z decay products Evidently, for this approach the relevant coupling is cj)ZZ 
(or (ftZ'-f). In this work we will consider a higher energy of 1 TeV for the e + e~ collider 
and examine the e~e + — > e~e + (f> process. At this energy VV fusion (V = Z,j) becomes 
dominant. Eventually, we will focus on the 77 fusion sub-process which gives a substantial 
production rate and could be analyzed to study the CP nature of the Higgs boson 0. 

In order to make any predictions for those processes involving the effective <pVV coupling, 
it is necessary to know the order of magnitude of the effective coupling constants g^yy. 
Several processes have been discussed in the literature that can be used to obtain bounds 
on these coefficients. In particular, based on the ELA, the constraints from the LEP and 
the Tevatron on the reactions e + e~ — > 777, pp — > jjjj, pp — > 77 + E T , and pp — > 777 
are known to yield bounds on fi/A 2 0, To be consistent with the bounds reported in 
Ref. {2I we will consider \fi\/A 2 < 50 TeV -2 , which implies |<70vv| < 0.3 in our notation, cf. 
Eq. JU). 

1 In general, the CP-odd Higgs boson can have a mass different from the CP-even Higgs boson. Since the 
goal of this study is to develop methods for determining the CP property of the Higgs boson, we assume 
a similar mass for these Higgs bosons to simplify our comparison. 
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III. MONTE CARLO ANALYSIS FOR e e+ -» e~e+0 



We now turn to discuss how the e~e + — > e~e + <fi process can be used to determine the 
CP property of the Higgs boson. In the following, H stands for the SM Higgs boson, and as 
discussed above, h and A denote the CP-even and the CP-odd Higgs boson, respectively. 

The Feynman diagrams for e~e + — > e~e + (p are shown in Fig. [TJ For the SM Higgs boson 
if, the contributions arise from Higgstrahlung (the left diagram) and vector-boson fusion (the 
right diagram) processes. For the CP-even Higgs boson h, in addition to the SM diagrams 
there would be contributions arising from the higher dimension hVV couplings. In contrast, 
for the CP-odd Higgs boson A, the contributions arise only from the higher dimension AVV 
operators, as there is no SM-like tree level AZZ coupling. 

We will carry out a Monte Carlo study and construct some CP-odd observables which are 
sensitive to the strength of the anomalous <pVV couplings and can be used to discriminate 
a CP-even from a CP-odd Higgs boson. In our analysis we will consider both unpolarized 
and polarized lepton beams, and the convention for the momenta of the external particles is 
taken as follows 

e~(pi)e + (p 2 ) -> e"(p 3 )e + (p 4 )0(p5) , (3) 

where p\ and p 2 are the incoming momenta and pz, p^ and p§ are out-going momenta. 

Given the effective interaction of Eq. (JTJ, it is straightforward to calculate the helicity 
amplitudes for e~e + — > e~e + 0. To this end we used the spinor product formalism introduced 
in Refs. Q Q- The resulting amplitudes are presented in Appendix A for the sake of 
completeness. To numerically illustrate the method of our analysis, we will take the Higgs 
boson mass to be 115 GeV (for each (f) — H,h, A) and the center-of-mass energy of the e + e~ 
collider to be 1 TeV. 

We will discuss first the case of unpolarized beams. 



A. Initial considerations and basic cuts for the 77 fusion i-channel diagram 

The e~e + — > e~e + <fi processes can be divided into two sub-processes: the Higgstrahlung 
process (V(p production with subsequent decay of the vector boson V into the pair of charged 
leptons) and the vector-boson fusion process. The Higgstrahlung process is s-channel like 
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(a) 




Figure 1: Feynman diagrams for e e + — ► e e + (/>: (a) the effective theory where the heavy dot 
denotes the effective 4>VV coupling (V = 7, Z); (b) the Standard Model. 

while the vector-boson fusion process is t-channel like. We can separate them by examining 
the invariant mass of the outgoing e~e + pair. (Of course, for Z(p production one can include 
the other decay modes of the Z boson which are not considered in this study). In the 
Higgstrahlung process, the two final state charged leptons originate from the vector-boson 
decay. Therefore, if we require the invariant mass (M 34 = m e + e - ) to be within a 10 GeV 
range of the Z-boson mass, the cross section will come mostly from the s-channel diagram. 
Otherwise, the main contribution will come from the t-channel (ZZ fusion) diagram. 

On the other hand, both s- and t-channel diagrams may have a singular behavior when 
the 077 (or 4>Z^) coupling is considered. In particular, collinear and nearly on-shell photon 
divergences appear, and we must impose appropriate cuts to avoid these kinematical regions. 
In this work we will start by imposing the following conditions: 

I cos 6*3 1 , 1 cos #4 1 < 0.99 , 

Pts,Pt4 > 10 GeV, (4) 
M 34 > 5 GeV . 

Here, the angle 63 (84) is the polar angle of the outgoing electron (positron). The incoming 
electron is moving in the +z axis direction. Pt denotes the transverse momentum. The first 
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Table I: The unpolarized cross sections of e~e + — ► e~e + 4>, where (j> is either the SM Higgs boson H, 
the CP-even boson h or the CP-odd A. We have taken y/S = ITeV, = 115 GeV, gnw = 0.1 
and gAW = 0.1 for V = Z or 7. (Only one anomalous coupling is taken to be non-zero for each 
case.) We have imposed the basic cuts of Eq. (J3J) in order to avoid the divergent collinear behavior 
of the scattering processes arising from the ^77 couplings. 



process 


cx(fb) 

\M 3 4-m z \ < 10 GeV 


<7(fb) 

|M 34 -mz| > 10 GeV 


e~e + — ► e~e + H 


0.42 


9.0 


e~e+ -> e~e + h {ghzz) 


0.88 


8.6 


e~e+ -> e~e+/i (5h 77 ) 


0.44 


13.9 


e"e + -> e"e + ^ («Mzz) 


0.21 


0.17 


e"e + -> e"e+A (<7A 77 ) 


0.02 


4.22 



two cuts are required to avoid collinear divergences in the t-channel process, whereas the third 
cut is required to avoid the divergence when the s-channel photon goes on-shell when ignoring 
the electron mass in the helicity amplitude calculations. Furthermore, the requirement on 
the cosine of the polar angles | cos^l < 0.99, which corresponds to 172° > 9 > 8°, will ensure 
that the final state leptons fall within the detector coverage at the ILC. 

The cross section for e~e + —> e~e + (f> with unpolarized beams and a fixed value of the 
coupling coefficients is shown in Table HJ The SM prediction for the e~e + — > ZH — > e~e + H 
process is about 0.4 fb, much smaller than the 9 fb predicted for the t-channel ZZ fusion 
process at y/S = 1 TeV. The purpose of this work is to distinguish between the CP-even 
and CP-odd scalars, and we will concentrate on studying the t-channel process. Table |U also 
shows that the contribution from the ^77 coupling is potentially larger than the contribution 
from the anomalous (j)ZZ coupling with the same magnitude. This is because the collinear 
divergences mentioned above only appear in scattering processes with nearly on-shell photon. 

Apart from analyzing the usual kinematic distributions, we will also consider the following 
angular variables which will be shown later to be sensitive to the CP nature of Higgs boson: 

X e - e + = X = (j3 3 - j3 4 ) ■ pi = cos 6» 3 - cos 6> 4 , 

v _ v P3 • (pi X p 4 ) p 4 ■ (pi X p 3 ) 

\Px x Pi\ \pi x pa] 
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where the unit vectors pi = f>i/\pi\. Generally, the experimental analysis of these variables 
would require particle identification of the outgoing leptons. However, the absolute values 
of X e - e + and Y e - e + could also be used with the advantage of not requiring to determine the 
charge of the outgoing leptons, separating e + from e~ in the final state, to discriminate a 
CP-even from a CP-odd Higgs boson. 

B. Comparison of production rates from the ^77, 4>Z^ and (j)ZZ couplings in the 
t-channel process 

For the rest of this work we focus our attention on the t-channel process, so that we will 
require the e~e + invariant mass to be away from the Z mass value: 

\M U - m z \ > 10 GeV. (6) 

It is understood that the conditions of Eqs. PJ and (JBJ) are being applied to the following 
analysis. We will refer to these conditions as basic cuts. 

In Fig. [2] we show the e~e + — > e~e + h production cross section induced by each of the 
3 couplings, (7h 77 , ghz-y and gtzz- Only one coupling is taken to be non-zero for each case. 
Because of the collinear divergent behavior originated from the photon propagator, the con- 
tribution from the ^77 coupling is much larger for similar values of the coupling coefficients. 

For the case of e"e + — > e~e + A production, the cross section depends directly on the gAvv 
coupling. We can write down the numerical value of the production cross section as: 

a(A 77 ) = 421 x f All (fb) , 

a(AZ^ = 145 x ~gl /Z (fb) , (7) 
a(AZZ) = 21Axf AZZ (fb) . 

In general, it is not possible to predict the relative sizes of the anomalous coupling coeffi- 
cients without knowing the underlying theory. In case that the underlying physics induce all 
three effective couplings with similar size, the ^77 coupling will produce the biggest effect 
on the t-channel production rate of Higgs boson. Thus, hereafter, we will assume that the 
effect of this coupling indeed dominates the Higgs boson production rate and we shall focus 
our attention on the effect of this coupling to various kinematic distributions of the CP-even 
and CP-odd Higgs boson signal events. 
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Figure 2: The unpolarized e~e + — ► e~e + h production cross section for different values of the 3 
CP-even coupling coefficients: <?h 77 , guz^ and ghzz- The basic cuts, Eqs. (J3J and (j^J, have been 
applied, and only one anomalous coupling is taken to be non-zero for each case. 

C. Unpolarized beams 

After imposing the basic cuts, cf. Eqs. (J4J and (JHJ), we show in Fig. 03 a few kinematic 
distributions which are useful for determining the CP property of the Higgs boson. They 
are the distributions of the transverse momentum and energy of the Higgs boson and lepton, 
and the invariant mass of the two final state leptons. 

The results shown in Fig. El can be understood as follows. The small tail in the low mass 
region of the lepton pair invariant mass plot suggests that the production rate is dominated 
by the ZZ and/or 77 fusion processes, and the s-channel process is not as important as 
the t-channel process in high energy collision with center-of-mass energy around 1 TeV. 
The SM Higgs boson production is through the dimension-4 HZZ interaction only. While 
the CP-even Higgs boson h could be produced through either the dimension-4 hZZ or the 
anomalous ^77 interaction, the CP-odd Higgs boson could only be produced through the 
A77 interaction. Because of the propagator of the Z-boson in the ZZ fusion process, the 
typical transverse momentum of the final state lepton after emitting the Z boson in the 
ZZ fusion process is about half of the Z boson mass. In case of the production of the 
CP-odd Higgs boson via 77 fusion, the final state lepton typically has a smaller transverse 
momentum, cf. Fig. Etc), due to the collinear enhancement in the emission of an effective 
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Figure 3: The normalized differential cross sections for the e~e + — > e~e + 4> processes as a function 
of the corresponding kinematical variable. We show the distributions for the SM H (blue dotted 
line), the CP-even h (black solid line) and the CP-odd A (red dashed line) Higgs bosons. The values 
g h ^ = c/A'yy = 0.1 have been used. The basic cuts ( Eqs. Q and (0) ) have been applied. 

photon from the lepton. Since the effective photon is transversely polarized, the lepton 
after emitting the effective photon in the 77 fusion process can have a large angle relative 
to the beam axis, and yields a smaller magnitude in the absolute value of its rapidity. In 
contrast, after emitting a longitudinal (virtual) Z boson, the lepton prefers to move along the 
forward direction which is obvious from the consideration of helicity conservation. The net 
effect of the collinear enhancement and the transverse polarization of the emitted effective 
photon to the lepton rapidity distribution can be seen in Fig. Ufa), which shows the rapidity 
distribution of the final state electron with the choice of the incoming electron in the +z 
direction. The distribution of the final state positron is not shown because it has the same 
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shape as the electron's but in the negative rapidity direction. As compared to the lepton 
rapidity distribution predicted by the SM, the one predicted for the CP-odd Higgs boson 
shows a flatter distribution. The result for the CP-even non-SM Higgs boson h lies between 
the above two curves. In Fig. Efb), we also show the distribution of the rapidity difference, 
in its absolute value, between the two final state leptons in the Higgs boson events. By 
taking only the absolute value of the rapidity difference between the two charged leptons, 
we eliminate the requirement of identifying the electrical charge of the two outgoing leptons, 
i.e., separating e + from e~, which could be challenging for leptons with large rapidity (in 
magnitude). For completeness, we also show in Fig. HJc) the distribution of the lepton 
rapidity difference without taking its absolute values (for the case that the electrical charge 
of the lepton is identified). Due to the collinear enhancement discussed above, very little 
event rate falls into the negative region. 

Next, let us discuss the kinematic distributions of the Higgs bosons. As shown in Figs.[SJa) 
and (b), the energy of the CP-odd Higgs boson is typically higher than that of the SM Higgs 
boson because the invariant mass of the final state lepton pairs in the CP-odd Higgs signal 
events is lower than that in the SM Higgs signal events, cf. Fig. Efd), which is due to the 
lower transverse momentum that the final state leptons can have in the CP-odd Higgs events. 
Furthermore, due to the transversity of the emitted effective photon in the 77 fusion process, 
the possibility for producing two energetic leptons from 77 fusion process is small, therefore 
the energy distribution of the CP-odd Higgs boson drops very fast in the small energy region, 
cf. Fig.^b). The same reasoning also explains the transverse momentum distribution of the 
Higgs boson in Fig.[Bta). In all cases, the maximal allowed value for the energy of the Higgs 
boson is constrained by the center-of-mass energy (\/S) of the collider and the masses of 
final state particles, it is 



y/S (m e ++m e -) 2 - 



711 



2 



EL* = V - — ^7= " ~ 507 GeV, (8) 

as can be deduced from the three-body phase space. 

We note that in the above kinematic distributions, the CP-odd Higgs boson differs from 
the CP-even Higgs boson, but the CP-even Higgs boson has similar distributions as the SM 
Higgs boson. To separate the CP-even Higgs boson from the SM Higgs boson background, 
one could make use of the rapidity distributions of the final state leptons, as shown in Fig. SI 
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Figure 4: The normalized differential cross section for the e~e + — ► e~e + (j) process as a function of 
the rapidity (difference). We show the distributions for the SM H (blue dotted line), the CP-even h 
(black solid line) and the CP-odd A (red dashed line) Higgs b osons. The values (7/177 — 9A-y^ — 

0.1 

have been used, and the basic cuts have been applied. 



Table II: The unpolarized cross sections of e~e + — > e~e + (p, where eft is either the SM Higgs boson 
H, the CP-even boson h or the CP-odd A. We have taken = 115 GeV and considered the 
following = ±0.1 and 5U77 — 0.1. The kinematic cuts listed in each column are applied 

sequentially. 



process 


(fb) 
basic cuts 


(fb) 

\Ve~ ~ Ve+ \ < 3 GeV 


(fb) 

\Ve--Ve+\ <2 GeV 


e~e + — > e~e + H 


9.0 


1.09 


0.17 


e~e + — > e~e + h (g^ 77 = 0.1) 


13.9 


3.50 


1.42 


e~e + — > e~e + h (5/177 = — 0.1) 


13.7 


3.30 


1.28 


e~e + — > e~e + A (^ 77 = 0.1) 


4.22 


1.91 


0.91 



Since the CP-even Higgs boson h can be produced either from ZZ or 77 fusion process, the 
electron rapidity distribution lies in between the SM Higgs boson and CP-odd Higgs boson. 
Based on the different shapes in the rapidity distributions of H, h and A, cf. Fig. HJ we 
can impose further cuts to reduce the SM rate (which is referred below as the "background 
rate"). In Table ITT1 we show the cross sections after demanding the rapidity difference between 
the two leptons to be smaller than 3 and 2, separately. The former cut, \q e - — r] e +\ < 3, 
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increases the signal-to-background ratio by about a factor 2 while keeping about 25% of the 
signal rate for the CP-even Higgs boson for both positive and negative g/i 77 . For the CP-odd 
Higgs boson, this cut increases the signal-to-background ratio by about a factor of 3.7 while 
keeping about 45% of the signal rate. That more CP-odd Higgs boson events survive this 
cut is due to the different electron rapidity distributions as discussed above. As shown in 
Fig. Ilfb), it would be possible to make a substantial reduction (by an order of magnitude) 
of the SM rate by imposing a lower rapidity difference cut. We observe, nevertheless, that 
for the stringent cut (\rj e - — r] e + \ < 2) with couplings and g^ 77 of order 0.1, only a small 
cross section of about less than lfb may be left. This may make it difficult to observe any 
events in a more realistic analysis. We therefore choose not to use this stronger cut. For the 
remainder of this work, in addition to the basic cuts we will apply the following cut 2 : 

\Ve- ~Ve+ \ < 3. (9) 

But it is worth mentioning that with a high luminosity in the upgraded ILC, one could use 
a tight cut to achieve a better significance. 

In Fig. we present a few kinematic distributions after imposing the lepton rapidity 
difference cut, cf. Eq. (0), together with the basic cuts. First of all, we note that the typical 
invariant mass of the lepton pair becomes smaller because the two leptons have a smaller 
rapidity separation. Furthermore, the SM Higgs boson rate is suppressed enough that the 
CP-odd Higgs boson rate becomes comparable to the SM rate. From the conservation of 
total energy, a smaller m e - e + implies a larger energy E$ (and transverse momentum p^,) of 
the Higgs boson which is evident in Fig. O Again, the m e - e + distribution of the CP-even 
Higgs boson h is similar to the CP-odd Higgs boson A in the small invariant mass m e - e + 
region while similar to the SM Higgs boson H in the large invariant mass m e - e + region. 
However, the detailed distribution depends on the interference pattern generated by the 
inclusion of the dimension-6 /177 operator. Clearly, the difference in the above mentioned 
kinematical distributions of the lepton pair and Higgs boson originates from different spin 
correlations in the scattering amplitudes for producing the three different Higgs bosons. To 
further investigate their difference, we also show in Fig. Elthe distributions of cos(0 e - — 9 e +) 

2 Notice that the cuts imposed in this work do not require identifying the electric charge of the leptons, i.e., 
separating e + from e~. 
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Figure 5: The normalized differential cross sections for the unpolarized e~e + — ► e~e + 4> processes as 
a function of the corresponding kinematical variable. We show the distributions for the SM H (blue 
dotted line), the CP-even h (black solid line) and the CP-odd A (red dashed line) Higgs bosons. 
The values = 5A77 =0.1 have been used. The basic cuts and the lepton rapidity difference 
cut, cf. Eq. (jHl, have been applied. 

and cos(0 e - — 4> e +), where Q{ and <f>i are the polar and azimuthal angles of the lepton % in 
the final state, and the two new observables X e ~ e + and Y e ~ e +, cf. Eq. The cause of the 
difference in the distributions of the above physical observables can be better understood 
when we study the case of polarized beams which will be discussed in the next section. 

Since the polarization of the incoming electron and positron beams would be possible at 
the ILC, it is worth analyzing this possibility. Below we will show that polarized beams are 
indeed useful for distinguishing a CP-even from a CP-odd Higgs boson. 
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Table III: Cross section for e^ 2 — > e~e + <p when the initial state leptons are polarized. L and 
R stand for a left-handed or right-handed polarized lepton beam. We have taken = 115 GeV, 
9/177 = 0.1 and g/^ = 0.1. 





e Ai e A 2 ~~ y e e+ H 


e Ai e A 2 ~~ * e e+ h 


e Ai e A 2 ~~ y e 




LL 


LR 


RL 


RR 


LL 


LR 


RL 


RR 


LL 


LR 


RL 


RR 


basic cuts 


8.65 


13.17 


5.54 


8.65 


17.16 


13.85 


7.47 


17.16 


4.66 


3.77 


3.77 


4.66 


with \rj e - — r] e + < 3 


1.06 


1.57 


0.67 


1.06 


4.93 


2.33 


1.85 


4.93 


2.21 


1.61 


1.61 


2.21 



D. Polarized beams 

Let us analyze to what extent the polarization of the collider beams could help us to 
discriminate between CP eigenstates of Higgs bosons. For simplicity, we assume a 100% 
degree of polarization. Needless to say that in practice the beam polarization will be only 
partially longitudinal or partially transverse. Using the rapidity cuts given in Eq. (JUJ) we 
calculated the rates for the eT^ e* — > e~e + 4> (with Ai,A 2 = LorR) processes, which are 
shown in Table Mil for all the possible polarization states of the initial leptons. The rates 
shown in the first row are with the basic cuts only, and the rates in the second row are with 
the additional cut on the lepton rapidity difference, cf. Eq. (JHJ). Since we assume CP is 
conserved in this study, the production rate of e^e^ is the same as the e^e\ rate. 

For the sake of illustration we will present below the kinematic distributions for left-left 
and left-right handed polarized beams, i.e. e^ej — > e~e + (f> and e^e^ — > e~e + 0, denoted as 
LL and LR, respectively. 

1. Left-left handed polarized beams 

Fig. |3 shows various kinematic distributions for e^e\ — * e~e + (f> when left-left handed 
polarized beams are used, in which both the basic cuts and the lepton rapidity difference 
cut, cf. Eq. Q, have been imposed. It is evident from the figure that the use of polarized 
beams enhances the differences between some distributions, like cos (9 e - —6 e +) and Y e - e + 
distributions, cf. Eq. We note that for this polarization state of the collider beams 
there is no s-channel contribution to the production process. The SM Higgs boson H is 
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Figure 6: Same as Fig. [HI for left-left polarized beams. 

produced from ZZ fusion and the CP-odd Higgs boson A is produced from 77 fusion, while 
the CP-even Higgs boson h can be produced via both ZZ and 77 fusion processes. Through 
a simple algebra, one can show that 

Y e - e + oc sin0 e - sin6> e + [cos (<p e - + <p e +) - cos (0 e - - <p e -)} , 

therefore, Y e - e + and cos (<p e - — <f) e +) provide almost the same information. To investigate 
the different shapes in the cos(0 e - — <p e +) and Y e - e + distributions between h and H pro- 
ductions, we show in Fig. the distributions of cos (6 e - — 6 e +) and Y e - e + after switching 
off the hZZ interaction in the h production, namely, only the /177 interaction is included 
in this comparison. It is clear that due to the different CP property of h and A bosons, 
the distributions of Y e - e + are very different for h and A bosons, cf. Fig. Efb), while the 
distribution of cos (6> e - — 9 e +) is almost the same for h and A bosons, cf. Fig. 0a). Hence, 
we conclude that the azimuthal angle distributions are more sensitive to the CP property of 
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Figure 7: Distributions of cos {9 e - — 6 e +) and Y e - e + after switching off the hZZ interaction in the 
h production, namely, only the hjj interaction is included in this comparison. 

the Higgs bosons. Particularly, in the Y e - e + distribution, a dip occurs around Y = for the 
CP-odd A boson production. In contrast, the rate is enhanced at Y = for the CP-even 
h boson production. Needless to say that after turning back on the hZZ interaction, the 
distributions of the h boson production are also affected by the interference effect between 
the hZZ and /177 interactions which results in the actual distributions shown in Fig. |B1 
Similar conclusion also holds for the cos(0 e - — e +) distributions. 

2. Left-right handed polarized beams 

Next, we consider the case that the electron beam is left-handedly polarized while the 
positron beam is right-handedly polarized. For this polarization state of the collider beams, 
both s- and t-channel processes contribute to the production of Higgs boson. In Fig. |H1 we 
show the similar kinematic distributions discussed above for comparison. One distinct feature 
in this set of distributions is that there are enhancement rates (spikes) when the energy 
and transverse momentum of the Higgs boson are near 500 GeV, cos(# e - — 9 e +) and 
cos(</> e - — 4> e +) are near 1, and both X e - e +and Y e - e + are near 0. All of these spikes originate 
from the s-channel process. In the s-channel production process, the invariant mass of the 
di-lepton pair can take a small value, i.e., the di-lepton pair is produced from the decay 
of Z boson or a low invariant mass of virtual photon. Obviously, the enhancement factor 
is larger for a di-lepton pair produced from a low mass virtual photon than an on-shell Z 
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Figure 8: Same as Fig. but for left-right polarized beams. 



boson. That explains why, for example, in the Y e - e + distribution, both h and A production 
show a strong peak than the H production in which only the dimension-4 HZZ coupling is 
considered. 



IV. CONCLUSIONS 

We study the production of Higgs boson (0) via the e~e + — > e~e + 4> process at a future 
International Linear Collider (ILC) withv^5 = 1 TeV. At this energy the Standard Model 
ZZ-fusion rate dominates over the associated ZH production rate. Here, we consider a 
class of theory models based on the effective (f>VV (with V = Z, 7) couplings generated by 
dimension-6 operators. We also show how to use the kinematical distributions of the final 
state leptons to discriminate the CP-even and CP-odd Higgs bosons produced at the ILC, 
when the ^77 coupling dominates the Higgs boson production rate. Assuming the Higgs 
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boson is identified, we have found that the kinematical distributions (and correlations) of 
the final state e~e + pair can be used to test the CP nature of the Higgs boson. In particular, 
we have found that requiring the absolute value of their rapidity difference to be within 3 can 
greatly reduce the SM rate so that the contributions from anomalous couplings are relatively 
enhanced and it becomes easier to use 077 couplings to determine the CP property of Higgs 
boson. One advantage of these cuts is that it is not necessary to determine the electric 
charge of the final state e~e + pair. (Particle identification could be challenging for leptons 
with large rapidity values (in magnitude).) Furthermore, we define two kinematical variables 
X e - e + and Y e - e +. The shapes of these distributions for a CP-even Higgs boson are different 
from those for a CP-odd Higgs boson. To distinguish a CP-odd Higgs boson from a CP-even 
Higgs boson produced via e~e + — > e~e + <fi process at a 1 TeV ILC, it is most effective to use a 
left-handedly polarized electron beam colliding a left-handedly polarized positron beam, and 
the observable with the biggest discriminating power is the Y e - e + distribution, as defined in 
Eq. ©. 
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Appendix A: HELICITY AMPLITUDES FOR e e+ e e^ 



„ :he calculation of the helicity amplitudes for e e + — > e e + 0. 
We use the method of Ref. 0, which breaks down the four-dimensional Dirac algebra 



In this appendix we outline t 
e use the method of Ref. 0, , 
into an equivalent two-dimensional one. In the Weyl basis, Dirac spinors have the form 

(Al) 




for fermions and 

V>± = < (A= _ } (A3) 



for antifermions. cj± = ^JE ± |p| and Xa/2 stands for the eigenvectors of the helicity operator 
p ■ a with eigenvalue A = 1 for spin-up fermions and A = — 1 for spin-down fermions: 

/ cos 0/2 \ f -e^sm9/2\ 

Xi/2 = I ' , X-i/2 =[ ' \. (A4) 

I sin 0/2 / V cos 9/2 I 

We will use the shorthand notation \pii> for x±i/2- In this basis, the Dirac matrices have 
the form 

( I (" (' "V (A5) 

\10) \aj ) \0 -lj 

where a,j are the 2x2 Pauli matrices. We can write 

(A6) 

where 7± = (l,=Fa). Let us denote the helicity amplitudes by A4(A e -, A e +, A e -, A e +), where 
Aj represents the helicity of the particle i. The convention for the momenta of the exter- 
nal particles is given in Eq. (jSJ). After a straightforward calculation, the following helicity 
amplitudes for e~e + —> e~e + 4> are obtained. 
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p + a ■ p l 


l-l 


' 


-a-p j 




l jL J 



1. SM Higgs boson 



The SM Higgs boson can be produced through the Higgstrahlung or ZZ-fusion pro- 
cess, as shown in Fig. The former is s-channel like but the latter is t-channel like. 
The helicity amplitudes Mf^ z (A e -, A e +, A e -, A e +) for the SM Higgs boson production via 



A^5M (+ + ++) = C fac - 2 ^ I « (i-| 2 +)(4+|3-) , (A7) 



M^ z ( ) = C /ac7 -5-=^^r« (i_|4+) (2+|3-) , (A8) 



M%& z {- + +-) = C faC / „ ^"T,T H -or (2-[4+) (3+ [I-) , (A9) 



M^f(+ - -+) = C fac — 2 '^;y K (2+|4-) (3-|l+) , (A10) 



V 

e 3 


)e + (p 4 ,A 4 )i7(p 5 ) 






/ 9 


9 \ / 9 9 \ 

-mz)\p\- m z) 






(Pi 


- m|)(p2 - m |) 




2 9Hzze z L e z R 


(P? 


-m 2 z)(P 2 m-m 2 z) 




2g HZZ e z e z R 




-m 2 z)(P 2 m- m z) 


/ 


2 gnzzefef. 



M^ z (- + -+) = C fac { *» u ^™ (4+|3-)(l-l2+) 



2g HZZ e z R e z 



(Pi -m 2 z )(p 2 m -m 



M^ Z (+-+-) = C iac \ ( 2 29H *ff R (l + l2-)<4-|3+) 

1 (P%-™z)(P 2 a- m z) 



(2-|3+> (4+|l-)>, 



(2-|3+)(4+|l-) S. (A12) 



(p 2 -m z ){p 2 m -m%) 



Here 



e R = —sw, 



and the SM H-Z-Z coupling is 
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PI P3 




in which sw = sin 6^ and c\y = cos 9w with 6w being the weak mixing angle. The common 
factor Cf ac is defined as 

and the momenta of the intermediate state particles are defined as p q = p 2 — p±, Pk = Pi — P3, 
Pi=Pi+ P2, and p m =P3 + Pa- 

2. CP-even Higgs Boson 

The CP-even Higgs boson gets both the SM-like contribution and the anomalous contri- 
bution, 

M even = MfJ + MhVV - 
v =1 ,z 

The SM-like contribution can easily obtained from Sec. lA II bv replacing the SM Higgs boson 
(H) with the CP-even Higgs boson (h). The anomalous contributions are listed as follows: 



M hVV (+ + ++) = C f J hVV 2 e Y R { 2( Pk ■ p q ) (1 + |4-) (2-|3+) - (l+l y g _|3+) (2-1 j/ fe+ !4-) 



PIP, 



(A13) 

'e v f 

^{ 2(p k -p a ) (1-|4+) (2 + |3-) - (2-\j/ a Al-) (3+| ^_|4+ 

PkPa 



M hVV ( ) = C f J hV \ e Y R {2{p k .p q ) (1-|4+) (2 + |3-) - (2-1 j/ q+ \l-) (3+| A„|4+) 



(A14) 
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■yV v 

M hVV (- + +-) = C s J hVV *\ 6r { 2( Pl ■ p m ) (2-|4+) (3 + |l-) - (2-| j/ m+ \l-) <3+| fi_\4+) \, 



PiP 



in 



(A15) 



M hVV (+ - -+) = C s J hVV ff R \ 2( Pl . Pm ) (2+|4-) (3-|l+) - (2+|^ m _|l+) <3-|^ + |4-> }, 

Pi Pm 



(A16) 

^h(p k ■ Pa ) (4+|3— > (1-I2+) - (2-| pUl-) (3-| ?/^|4- 



+ = C /ac ^f^ <| 2{p k .p q ) (4+|3-> (1-|2+) - (2-UJ1-) (3-|^ + [4-) 



(A17) 

+ ^ 2 fo ' Pm) (2 " |3+) (4+|1 " ) " < 2 "l ^l 1 ") (3-1 ^ + l 4 -> \> 

Pi Pn 



(A18) 



M hVV (+-+-) = C f J hVV f L f R { 2( Pk ■ p q ) (1 + |2-) (4-|3+) - * g j3+) <2+| A_|4+) 

+ Cj J hyv f^ R \ 2(p, ■ p m ) (2 + |3-) (4-|l+) - <2+| ^ m _|l+> (3+| y,_|4+) j> . 



Here, 

1 , 1 



, + 40 v = z 

Cvi/ z 

-% V = 7 

^ \/ = Z 

-Sly V = 7 



and the anomalous couplings Qhvv are related to the coefficients /j in Eq. ((2J) by 

2 

Sy/ m W , t f \ 

gftVV — — g \JBB + /WW] • 

Again, we have defined p q = p 2 - p 4 , Pfc = Pi ~ P3, Pi = Pi + P2, and p m =Pi+p 3 - 



3. CP-odd Higgs Boson 



The CP-odd Higgs boson receives contributions only from the anomalous couplings, i.e. 
M odd = £ A4 A ^. 
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Using the method of Ref. |3j, we obtained the following helicity amplitudes: 



M AVV (+ + ++) = C, J AVV 2 e f R I (1 + |4-) <2-| A + AJ3+) - (3-|2+) (4+| A- A + |l-> 



3 



A^ VV ( ) = g/ac ^gP { (3+|2-) (4-| y >+ ^Jl+) - (1-|4+) (2+| A. ^ + |3-) 



M AVV '(- + +-) = C f J A ™ e } e * {(l+\3-) (4-1 y l+ ^ m _[2+) - (2-[4+) (3+1 jrf t _ 



= C S J AVV *. Y R \ (2+14-) (3-1 y l+ JC-11+) ~ (1-13+) (4f[ y m+ !2~) 



-M^(- + -+) = C/ J^rf ^ <3+[4-) (2-| ^ + ^Jl+) - (1-12+) (4+| $ k _ y g+ [3-) 



Pit? 



+ C/oc ^f 6 ^ < 1 + l 4 -) < 3 "l ^-|2+> - (2-|3+) <4+| jW|1-> 



M AVV '(+-+-) = C/a /^TP^ (1+12-) (4-1 A + A.I3+) - (3-14+) (2+1 

+ C s J AVy 2 e Y R { (2+|3-) <4-| # + ?f m _|l+) - (1-|4+) <3+| j^ + |2-) 
where the anomalous couplings cjavv are related to the coefficients fi in Eq. (j2J) by 

,2 



<Mvv — — 9 ^ 2 yjBB + /ww J ■ 
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